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ABSTRACT 

The advection dominated accretion flow (ADAF) has been quite successful in ex- 
plaining a wide variety of accretion-powered astronomical sources. The physical char- 
acteristics of ADAF complement the classical thin disk flow quite nicely, and extensive 
work has been done on it. However, all high temperature accretion solutions includ- 
ing ADAF are physically thick, so outgoing radiation interacts with the incoming flow. 
Thus, ADAF solutions share as much or more resemblance with classical spherical ac- 
cretion flows as with disk flows. This interaction, which has been neglected by most 
authors, is primarily through the Compton heating process which will typically limit 
the steady solutions to L < 10~^Le, where Le is the Eddington luminosity. We ex- 
amine this interaction for the popular ADAF (Narayan & Yi 1995a) case but similar 
conclusions, we expect, would apply for other high temperature, geometrically thick 
disks as well. We study the global thermal nature of the flow, with special considera- 
tion given to various cooling and, especially, preheating by Comptonizing hot photons 
produced at smaller radii. We flnd that without allowance for Compton preheating, a 
very restricted domain of ADAF solution is permitted and with Compton preheating 
included a new high temperature PADAF branch appears in the solution space. In the 
absence of preheating, high temperature flows do not exist when the mass accretion 
rate m = Mc^/Le ^ 10^^'^. Below this mass accretion rate, a roughly conical region 
around the hole cannot sustain high temperature ions and electrons for all flows having 
m > 10"^, which may lead to a funnel possibly filled with a tenuous hot outgoing wind. 
If the flow starts at large radii with the usual equilibrium temperature ~ 10^ K, the crit- 
ical mass accretion rate is much lower, m ~ 10~^'^ above which level no self-consistent 
ADAF (without preheating) can exist. However, above this critical mass accretion rate, 
the flow can be self-consistently maintained at high temperature if Compton preheating 



- 2 - 



is considered. These solutions constitute a new branch of solutions as in spherical accre- 
tion flows. High temperature PADAF flows can exist above the critical mass accretion 
rate in addition to the usual cold thin disk solutions. Therefore, the Compton preheat- 
ing could be the mechanism to trigger the phase change from the thin disk to PADAF. 
We also find solutions where the flow near the equatorial plane accretes normally while 
the flow near the pole is overheated by Compton preheating, possibly becoming, a polar 
wind, solutions which we designate WADAF. 

Subject headings: accretion, black hole, X-ray sources, QSOs 



Gas accretes onto compact objects in disk or (quasi) spherical shape depending on the angular 
momentum and the entropy of the gas. The flow may have a thin disk shape with very small 
radial motion, or it could be spherical or spheroidal with significant radial motion. The amount 
of radiation emitted depends on the shape of the flow as well as on the nature of the compact 
object. Accretion onto white dwarfs and neutron stars always produces luminosity approximately 
proportional to the mass accretion rate. Thin disk accretion also produces radiation that is a 
constant fraction of the gravitational energy. However, spheroidal or spherical flow onto the black 
holes with relatively large radial velocities emits differing amount of radiation, roughly scaling with 
the square of the accretion rate (Shapiro 1973a; Park 1990a,b), depending on the physical states 
of the gas and how the radiation escapes from the flow. Therefore, accretion onto black holes can 
have a variety of geometrical flow shapes and can produce a wide range of luminosity and emitted 
spectrum. 



Accretion flows onto black holes become spherical or nearly spherical when the angular mo- 
mentum is very small or when the hole is rapidly moving relative to the gas (Hoylc & Lyttlcton 
1939; Bondi & Hoyle 1944; Bondi 1952; Loeb &: Laor 1992). We review first the simpler case of 
spherical accretion because we find that the type of solution obtained here are quite relevant to 
more complex ADAF case. As long as two-body processes are important, the solutions become 
scale- free and are applicable to black holes of arbitrary mass (Chang & Ostriker 1985). These 
scale-free solutions then depend primarily on the dimensionless luminosity 



1. 



Introduction 



1.1. Spherical accretion 




(1) 
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and the dimensionless mass accretion rate^ 



M Mc^ 

Me Le ^ ^ 

or, equivalently, the dimensionless luminosity I and the radiation efficiency 

e=J^, (3) 

where is the Eddington luminosity (Ostriker et al. 1976). However, if we require the self- 
consistency between the gas and the radiation field generated by the accretion process, only certain 
combinations of the accretion rate and the luminosity (plus spectrum) are allowed, and a given 
self-consistent solution is described by a line in the {I, m) plane. Within a given type of solution, 
only one parameter is needed to characterize the specific flow, and the dimensionless mass accretion 
rate, rh, is the most meaningful parameter. Broad reviews of accretion physics can be found, for 
example, in Rybicki & Lightman (1979), Pringle (1981), Treves, Maraschi, & Abramowicz (1988), 
Prank, King, & Raine (1992), Chakrabarti (1996b). 



1.1.1. Spherical flow with small mass accretion rate (m <C Ij 

The electron scattering optical depth from infinity down to the horizon is roughly given by 
Tes = m{r/rs)~^^'^ for freely falling flow where is the Schwarzschild radius of the hole and, 
therefore, for m <C 1, spherical flow is optically thin to scattering. Photons escape through the flow 
almost freely, and the exchange of momentum and energy between the gas and the radiation field 
is negligible. Radiative cooling and Compton heating are also quite negligible. Gas is almost freely 
falling, v/c ~ (r/r^)"^/^, which is true even in relativistic flow (Shapiro 1973a,b, 1974); adiabatic 
heating keeps the gas up to the virial temperature, T oc (r/r^)"^, until electrons become relativistic 
(Shapiro 1973a; Park 1990a,b). Once electrons reach relativistic temperature, the adiabatic index 
of the gas changes from 7 = 5/3 to 7 = 13/9 if electrons and protons are well coupled. Now the 
temperature is increasing less rapidly and radiative cooling, mainly a relativistic bremsstrahlung 
and synchrotron, becomes more efficient. Regardless, the gas can reach as high as a few times 10^° K 
even when electrons are loosely coupled to ions via Coulomb process. The luminosity and efficiency 
of the non magnetic solutions are quite small due to the very low gas density, typically I < 10~^ 
and e < 10""^ {solid line marked as S [Shapiro 1973a] and large crosses [Park 1990b] in Figs. 1 
&; 3) with efficiency approximately proportional to the mass accretion rate e oc m and luminosity 
I oc m^. The radiation temperatures Tx, defined as 4:kTx being equal to the energy- weighted mean 
photon energy, of the solutions are plotted as crosses in Figure 2. 



^We define the Eddington mass accretion rate as Me = Le/c^ , without e ^ factor because the radiation efficiency 
e of accretion onto black holes is not fixed in general. 
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1.1.2. Spherical flow with large mass accretion rate (rh^ 1) 

As the mass accretion rate m approaches 1, radiative cooUng becomes more efficient due to 
the increased density and the flow cools down to ~ 10^ K in the absence of Compton heating. The 
flow forms an effectively optically thick, i.e., optically thick to absorption, core with blackbody 
radiation field inside. The flow no longer behaves as an adiabatic gas. 

Further increase of the mass accretion rate makes the core grow larger, and soon the whole 
flow becomes eflFectively optically thick. The flow is now cooled by the diffusion of radiation as 
in stellar interior (Flammang 1982, 1984; Soffel 1982; Blondin 1986; Park 1990a; Nobili, Turolla, 
& Zampieri 1991). However, the ultimate source of the luminosity is gravity through adiabatic 
compression rather than nuclear burning. Triangles in Figures 1-3 show the luminosity I, the 
radiation temperature Tx, and the radiation efficiency e as a function of the mass accretion rate 
m (Nobili et al. 1991). Note that the high temperature, low rh solutions {large crosses) are quite 
distinct from the low temperature, high rh solutions {triangles) in Figure 2, but they are continuous 
in Figure 1. Although the luminosity increases steadily with larger mass accretion rate, the flow is 
still a quite inefficient radiator, e < 10~^. 

The electron scattering optical depth becomes larger than unity when m > 1, and the flow 
becomes optically thick to scattering. Photons collide more often with the flow and the diffusion 
speed of photons decreases with increasing optical depth. For high rh, the bulk velocity of the flow 
exceeds the photon diffusion speed, v > c/{Tes + Tabs)-, and radiation trapping occurs (Begelman 
1978). The behavior of radiation and gas in this regime has to be treated under rigorous relativistic 
framework (Flammang 1982, 1984; Park 1990a; Nobih et al. 1991; Park 1993). In such an optically 
thick relativistic flow, the flux seen by the stationary observer relative to background spacetime is 
significantly different from the flux seen by the observer comoving with the flow. The momentum 
transferred to the gas from the radiation is closely linked to the former, a comoving-frame flux, 
and the luminosity seen by the observer at infinity to the latter, a fixed-frame flux. (Mihalas and 
Mihalas 1984; Park 1993). Also in this regime, the usual diffusion relation between the flux and 
the radiation pressure does not hold and should not be used (Miller 1990; Park &; Miller 1991). 
This radiation trapping critically affects all types of accretion flow with large mass accretion rates. 

The possibility of another mode at the same mass accretion rate was first investigated by Wan- 
del, Yahil, & Milgrom (1984). Subsequent rigorous studies show that the flow can be maintained 
up to high temperature through Compton preheating by the radiation that accretion flow itself 
produces at smaller radii (Park 1990a,b; Nobili et al. 1991; Mason &; Turolla 1992). They find this 
high-temperature, higher-luminosity, and therefore higher-efficiency flow does exist, but only for 
3 i$ 'Ti < 100 as shown in Figures 1-3 as large open circles (Park 1990b). Hence, for a certain 
range of the mass accretion rate, there exists more than one family of accretion flow with different 
properties. 

These high-temperature accretion flows have the gas temperature as high as 10^ ~ 10^^ K 
in contrast to the low-temperature (10'^ K) case. The flow is optically thick to scattering, and 
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radiation is trapped inside r ~ mr^. Yet, the flow is optically thin to absorption despite the 

large scattering optical depth because of the high electron temperature. The flows are much more 
luminous / ~ 10~^ — 10~^ {large open circles in Fig. 1) than their low-temperature counterparts 
(triangles) and produce much harder photons. The radiation efficiency e ~ 10"*^ is about 10^ times 
higher than for the low-temperature flow, yet still far less than e ~ 0.1 of the thin disk. The self- 
consistent luminosity I increases with the mass accretion rate m, while the spectrum of the emitted 
radiation which preheats the flow in turn is harder for lower rh and softer for higher rh {circles in 
Fig. 2). The global stabilities of these flows have been investigated in relativistic framework, and 
they are found to be thermally unstable and develop moving shocks (Zampieri, Miller, &: Turolla. 
1996) as suspected in steady-state work (Park 1990a,b). 

There exist a few candidate processes that might significantly increase the radiation efficiency. 
Meszaros (1975) considered the dissipational heating from turbulent motion, Maraschi et al. (1982) 
magnetic field reconnection, and Park &; Ostriker (1989) electron-positron pair production. Solu- 
tions of the first and second types can reach e as high as ~ 0.1, small squares in Figures 1 & 3 
(Maraschi et al. 1982), while that of the last type up to e ~ 10~^'^, small filled circles in Figures 
1 & 3 (Park & Ostriker 1989). They constitute yet another type of accretion flow in this so called 
super-Eddington accretion regime. However, the foundations for these higher-efficiency families are 
less solid than the other two families because the solutions are found without the consideration 
of preheating (Meszaros 1975; Maraschi et al. 1982) or only the inner supersonic part of the flow 
is dynamically considered (Park & Ostriker 1989). The rigorous and full treatment of dynamics 
and preheating will certainly show that they are unstable to the preheating instability (Cowie, 
Ostriker, & Stark 1978; Ciotti & Ostriker 1999). It is possible that the overheated flow is stabilized 
by developing a shock (Chang & Ostriker 1989; small crosses in Figs. 1 & 3), yet no self-consistent 
solutions with preheating and a steady shock have been constructed despite some hints of their 
existence (Park 1990b; Nobih et al. 1991). 



1.2. Preheating 

For this higher rh and I flow, radiation has more chance to interact with the gas: both momen- 
tum and energy are transferred from one to the other. Transfer of momentum, when the interaction 

of radiation with matter is mediated by the electron Thomson scattering, leads to the critical Ed- 
dington luminosity (/ = 1 lines in Figs. 1 & 3). Another lower type of luminosity upper limit is 
found by Ostriker et al. (1976) that can be as small as 10^^ times the Eddington luminosity. This 
new critical luminosity exists due to preheating of gas flow by hard Compton radiation. When the 
gas is preheated significantly near the sonic point, it becomes too hot to accrete in a steady state. 
For example, when the surrounding gas temperature is near 10^ K and the Compton temperature 
of preheating radiation is 10^ K, steady-state accretion is impossible above Icr — 0.01 {rectangular 
region I in Figs. 1 & 3). This preheating explains why self-consistent solutions are not found for 
m < 3 and rh > 100 (Park 1990a,b; Nobili et al. 1991): since Compton preheating depends on 
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the energy and the luminosity of outcoming photons, the preheating instabihty occurs when the 

radiation is too hard (m < 3; Fig. 2) or the luminosity is too high (m > 100; Fig. 1). The region 
can be smaller under certain conditions (Ostriker et al. 1976; Cowie et al. 1978; Bisnovatyi-Kogan 
& Bhnikov 1980; SteUingwcrf k Buff 1982; Krolik & London 1983). 

Preheated flow shows interesting time-dependent behavior (Cowie ct al. 1978; Grindlay 1978) 
that may be applicable to variabilities of accretion powered astronomical sources, including galactic 
X-ray sources, active galactic nuclei, QSOs, and cooling flows (Ciotti & Ostriker 1997, 1999). Cowie 
et al. (1978) found two distinct type of time-dependent behavior inside and near the sonic radius. 
The flow is overheated inside the sonic point for high I, low e accretion (region II in Figs. 1 & 
3), producing recurrent flaring on short time intervals. Low /, high e flow (region III in Figs. 1 & 
3) develops preheating outside the sonic radius, causing accretion rate and luminosity changes on 
longer time scales. However, these time-dependent studies did not impose self-consistency, meaning 
the radiation efficiency is arbitrarily assigned rather than calculated from the flow itself. 

All of the high temperature, physically thick self-consistent solutions shown in Figure 3 that 
are known or expected to be stable have a quite low efficiency (e < 10^^). It seems likely that all 
physically thick solutions having a high enough density to efficiently produce radiation will undergo 
relaxation oscillations with the primary cause being preheating instabilities. The exception to this 
may be the slim disks which have a relatively low emission temperature and which essentially 
resemble massive stars. 

1.3. Axisymmetric accretion 

Accretion flow onto a black hole with angular momentum can assume disk-like or spheroidal 
shape depending on how efficiently the heat produced is removed. The flow becomes thin disk- 
like when the energy from viscous dissipation is readily radiated away in the vertical direction 
(Shakura &i Sunyaev 1973). If too much radiation is produced and trapped within the flow, the 
disk puffs up to a thick one, and much of the radiation and gas energy is advected into the hole with 
significant radial velocity (Jaroszyiiski, Abramowicz, & Paczyhski 1980; Paczyhski &: Wiita 1980). 
When the disk is not thin, constructing the flow solution becomes a two-dimensional problem which 
has been tackled only by arbitrarily fixing certain physical parameters (e.g., Paczyhski 1998), or 
by reducing it to one-dimensional problem (Abramowicz et al. 1988; Narayan & Yi 1994, NYl 
hereafter; Narayan &; Yi 1995b, NY3 hereafter), or by numerical simulations (Molteni, Lanzafame, 
&; Chakrabarti 1994; Ryu et al. 1995; Chen et al. 1997; Igumenshchev & Bcloborodov 1997; 
Igumenshchev &: Abramowicz 1999; Stone, Pringle, &; Begelman 1999), or by assuming self-similar 
forms (Narayan Sz Yi 1995a, NY2 hereafter). 
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1.3.1. Thin disk 

Since the seminal works of 70's, thin disk accretion, especiahy the so-called a-disk has been de 
facto standard accretion mode applied to various accretion powered sources (Shakura 1972; Pringle 
& Rees 1972; Shakura & Sunyaev 1973; Novikov & Thorne 1973). The thin disk has merits of 
being simple: the equations are reduced to one-dimensional ones and all physical interactions can 
be described by local quantities. Angular momentum of the rotating gas is transported outward 
by the viscous stress through the differentially rotating gas. Gravitational potential energy is 
radiated away via viscous dissipation. Although the disk is geometrically thin, it is optically thick 
to absorption, and radiation diffuses out in the vertical direction. Because of the geometrical shape, 
the emitted radiation does not in general interact with other parts of the flow. The innermost edge 
of the disk determines the amount of gravitational potential energy to be released, and the radiation 
efficiency is roughly 0.1 regardless of the central compact object being neutron star or black hole. 
Therefore, all thin disk accretions appear on the dotted line e = 0.1 (TD in Figs. 1 &: 3). 

However, the thin disk formalism may not be extended to high luminosity systems or, equally, 
high m systems: the disk becomes unstable to thermal or secular instability when radiation pressure 
dominates over the gas pressure (Lightman &: Eardley 1974; Shakura & Sunyaev 1976; Pringle 1976; 
Piran 1978). Furthermore, hard X-ray sources like Cygnus X-1 cannot be explained by a thin disk: 
it is too cool and emits blackbody-like spectrum. Though Shapiro, Lightman, &; Eardley (1976) 
found a new family of optically thin disk with very high electron temperature, they are also proved 
to be unstable to the thermal instability (Pringle 1976; Piran 1978; Park 1995). Thus, for high flow 
rates or high temperature flows, this simple, high efficiency solution is not available or appropriate. 

1.3.2. Thick disk 

As the dimensionless mass accretion rate m approaches (i-e., ^ ^ 1), the vertical height 
of the disk becomes comparable to the radius. In the thin disk case, gravity is balanced by the 
Keplerian rotation and the radial motion is negligible. When the disk becomes geometrically thick, 
the dynamical balance in the radial direction becomes important, and significant radial motion 
occurs, resulting in sub-Keplerian rotation. Now the disk is two-dimensional: dynamics in the 
radial direction as well as in the azimuthal direction should be considered. 

This geometrically thick accretion disk has been studied by Paczyhski and collaborators (e.g. 
Paczyriski & Wiita 1980; Jaroszynski et al. 1980; Abramowicz, Calvani, &: Nobili 1980). Thick 
disk flow can only be described by full two-dimensional partial differential equations because the 
angular momentum of the gas is not known a priori, radial infall motion can not be ignored, and 
the energy equations become non-local. Studies of thick disk onto black holes show that the flow 
develops a cuspy inner edge through which the gas flows into the hole with roughly free-fall velocity. 
Photons are locked inside the disk and advected into the hole; yet a large amount of radiation may 
escape through the surface of the thick disk which shapes up the funnel around the rotation axis. 
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1.3.3. Slim disk 

To eschew the complexities while retaining the essential features of the thick accretion flow, 
Abramowicz et al. (1988) applied height-integrated a disk equations to super-Eddington (m > e~^) 
accretion flow. They found stable disk solutions with significant radial motion in high m regime, 
in which the thin disk is unstable. The flow is optically thick and the radiation is trapped at 
small radii so that, as noted above, these solutions resemble massive stars. The entropy is directly 
advected with the flow, which adds to the radiative cooling through the surface of the disk. This 
new family of disk solutions is called 'slim disk' because the disk is assumed to be slim enough to 
validate the height-integrated equations as in the thin disk solutions. In this accretion mode the 
radiation efficiency is not predetermined: some of the gravitational potential energy is ultimately 
absorbed into the hole due to the entropy advection, and the amount of radiation emitted depends 
on specific conditions of the flow. In this regard, the slim disk has some similarity to the case of 
high rh spherical accretion. The luminosity and radiation efficiency of these slim disks are plotted in 
Figures 1 & 3 as dot-dashed curve labeled SD; the flow can reach super-Eddington luminosity, I > 1, 
and e decreases from 0.1 as m increases (Szuszkiewicz, Malkan, & Abramowicz 1996). However, it 
is not at all clear if the slim disk approach can be extended to super-Eddington systems because 
the disk is likely to be thick rather than slim. 



1.3.4- Advection dominated accretion flow 

The slim disk solutions were constructed by explicitly integrating the height-averaged hydro- 
dynamic equations (Abramowicz et al. 1988). However, Narayan & Yi (NYl) find that these 
one-dimensional equations admit simple self-similar solutions if cooling due to the entropy advec- 
tion is assumed to be a constant fraction of the viscous heating (also see Spruit et al. 1987): The 
density, velocity, angular velocity, and total pressure are simple power laws in radius: p oc r~"^/^, 
V oc r~^^^, ft oc r^/^, and P oc r~^^^. Most of the dissipation energy generated is advected in- 
ward with the flow, and the solutions are called advection dominated accretion flow ^ (ADAF; see 
Narayan, Mahadeva, & Quatacrt 1998b for review). This type of accretion had been studied earlier 
(Ichimaru 1977) under different names like 'ion torus' (Rees et al. 1982). 

However, since the temperature of this type of flow is close to the virial value, it is not certain 
whether height-integrated equations are adequate. In the subsequent work (NY2), the full two- 
dimensional self-similar solutions have been constructed under the standard a- viscosity assumption. 
The radial velocity of the flow turns out to be a fraction of (roughly proportional to a for a ^ 1) 
the free-fall value at the equatorial plane, yet the flow along the pole is pressure supported with 
zero infall velocity. The matter is preferentially accreted along the equatorial plane. The flow is 



^The basic formulation of ADAF is the same as that of the slim disk flow. They are different only in m, or in the 
effective optical depth. 
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subsonic due to the near virial pressure, which is vahd for r 3> r^. Another unique feature of ADAF 
is that the Bernoulh constant of the flow is positive, as noted by NY2 themselves and critically 
analyzed by Blandford & Begelman (1999), which has led Xu & Chen (1997) and Blandford & 
Begelman (1999) to propose a generalization of ADAF that has both inflow and outflow at the 
same time. Another important result from NY2 is that the physical quantities calculated from the 
vertically integrated one-dimensional equations generally agree with those averaged over the polar 
angle. In a way this justiflcs the slim disk approach to the two-dimensional accretion flow, althoTigh 
this is proved only in the self-similar flow. This permitted further extensive works on the dynamics 
of ADAF to be based on the one-dimensional height-integrated hydrodynamic equations (Narayan, 
Kato, & Honma 1997; Chen, Abramowicz, & Lasota 1997; Nakamura et al. 1997; Gammie k, 
Popham 1998; Popham &; Gammie 1998; Chakrabarti 1996a; Lu, Gu, & Yuan 1999). Because, 
everywhere in the flow the radial velocities are quite substantial and the pole to equator density 
variations are moderate, these solutions share many of the characteristics of the purely spherical 
flow. 

Detailed microphysics are readily incorporated into these one-dimensional solutions, showing 
the main characteristics of ADAF: low radiation efficiency and high electron temperature (NY3; 
Abramowicz et al. 1995). These make ADAF applicable to low luminosity, hard X-ray sources, 
like Sgr A* (Narayan, Yi, & Mahadevan 1995; Mahadevan, Narayan, & Krolik 1997; Manmoto, 
Mineshige, & Kusunose 1997; Narayan et al. 1998a), NGC 4258 (Lasota et al. 1996), soft X- 
ray transient sources (Narayan, McClintock, & Yi 1996; Narayan, Barret, &; McClintock 1997), 
and low-luminosity galactic nuclei (Di Matteo & Fabian 1997; Mahadevan 1997) as well as to other 
diverse problems like lithium production in AGN (Yi & Narayan 1997), torque-reversal in accretion- 
powered X-ray pulsars (Yi & Wheeler 1998), the QSO evolution (Yi 1996; Choi, Yang, & Yi 1999), 
and X-ray background (Yi & Boughn 1998). Constructing the accretion flow with the inner part 
being ADAF and the outer part being thin disk has more merits and freedom than the simple thin 
disk, and therefore it can be successful in representing a wider range of astronomical systems. 

Typical ADAFs are presented as dashed lines in Figures 1 & 3 (NY3). They are inefficient 
radiators: most of the gravitational energy of the accreted matter is carried with the flow into the 
hole without being released in the form of radiation. This very property of ADAF is similar to 
that of the hot, adiabatic, low-m spherical accretion flow (Shapiro 1973a, b) as the slim disk is to 
the low-temperature, high-m spherical flow (Flammang 1982, 1984). The self-similar functions of 
radius for various physical quantities are exactly those of the non-relativistic adiabatic spherical 
flow, and the isodensity surface of the ADAF is spheroidal rather than disk-like (NY3). In fact, 
optically thin two-temperature spherical accretion flow with a magnetic field produces a spectrum 
that roughly agrees with that of ADAF (Melia 1992, 1994). 
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1.3.5. Unified description 

Recent works on various types of axisymmetric accretion flow with a wide range of mass 
accretion rate make it possible to describe these seemingly distinct families of accretion flow in a 
unifled way (Abramowicz et al. 1995; Chen et al. 1995). Each family of solution is constructed on 
appropriate assumptions, and specific name is attached to it, e.g., thin disk, slim disk, and ADAF. 
The traditional way of presenting the disk solutions in a plane of the mass accretion rate, m, versus 
the surface mass density, S, at a given radius is very revealing, contrast to the spherical accretion 
where the global properties like the mass accretion rate and the total luminosity are used. Different 
families of solutions appear as different curves in this plane (Fig. 4) . 

Since the solutions sometimes depends on the mass of the hole M and the value of the viscous 
parameter a, we take a specific case of M = 10 Af© and a = 0.1 (Chen et al. 1995). Figure 4 shows 
two disconnected curves, each representing the low-S flow {dotted curve on the left) and the high-S 
flow {solid curve on the right). The former is optically thin to absorption and the latter optically 
thick. This shows that there can be more than one type of accretion flow in a certain range of mass 
accretion rate. Also, the former type of accretion flow exists only for m <C 1, while the latter type 
is possible for all range of m. 

The lower (GTD) and middle (RTD) branches of the S shaped curve on the right side of Figure 
4 correspond respectively to the classic gas pressure and radiation pressure dominated thin disk 
solutions (Shakura & Sunyaev 1973). Radiative cooling through the surface of the disk is dominant 
over the advective transport. The positive slope of the gas pressure dominated disk shows the flow 
being stable, and vice versa for the radiation dominated disk (Lightman & Eardley 1974; Shakura 
& Sunyaev 1976). The uppermost branch (SD) of the curve is the slim disk solution in which the 
advection of radiation plus gas stabilizes the radiation pressure dominated flow. The vertical scale 
height of the slim disk can be comparable to the radius, especially when I = em — > 1. Even at 
this high m, the temperature of the flow is less than ~ 10^ K, and the spectrum of the emitted 
radiation is the superposition of modified blackbodies with different temperatures. Both the gas 
prcssTU'c dominated and the radiation pressure dominated thin disks are efficient radiators e ~ 0.1 
(TD in Figs. 1 &; 3) while the slim disk can have smaller efficiency e < 0.1 (SD in Figs. 1 & 3). 

The right branch (SLE in Fig. 4) of the optically thin accretion disk is the two-temperature 

hot disk solution of Shapiro et al. (1976). The slope indicates that the disk is stable to viscous 
perturbations. However, it is thermally unstable on much shorter time scales (Pringle 1976; Piran 
1978). The fiow has different ion and electron temperatures, and the viscous dissipation energy 
is readily radiated through the optically thin flow. Therefore, the radiation efficiency is always 
e ~ 0.1. Recent work on the stability of SLE solutions suggests that the advective cooling may 
affect the thermal instability (Wu 1997). 

The left branch (ADAF in Fig. 4) represents the ADAF of NYl and Abramowicz et al. (1995), 
in which the flow is optically thin and viscous dissipation energy is mostly advected with the flow, 
therefore, stabilizing the flow. The electron temperature of the flow can reach up to ~ 10^ K, and 
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the radiation spectrum is that of the Comptonized bremsstrahlung and synchrotron. The radiation 
efficiency of this branch (NY3) is roughly e oc m (see NY in Fig. 3) and the luminosity I oc (see 
NY in Fig. 1) as in the optically thin spherical accretion (NY3). 

There also can be other interesting families of axisymmetric accretion if processes like pair 
production (Kusunose & Mincshige 1992) or various shocks (Chakrabarti 1996a) exist, which have 
properties analogous to the equivalent solutions in the spherical flow. However, preheating, which 
is important in the spherical flow, has never been incorporated in any of the disk solutions shown 
in Figure 4. 

1.4. Two-dimensional nature of axisymmetric accretion flow 

One important property of the slim or thick axisymmetric accretion flows, including ADAFs, is 
their two-dimensional nature. This has been generally ignored under the assumption that height- 
integrated equations are good enough (NY2). However, Park & Ostriker (1999, hereafter PO) 
emphasize that the radiation emitted must interact with the accreting matter due to the geometrical 
shape of ADAF . Unlike the thin disk where the radiation easily escapes the flow through the vertical 
surface without any interactions, photons produced at smaller radii of the flow must escape through 
the tenuous outer part of the flow in much the same way as in spherical accretion. By studying 
the interaction of radiation with matter in two-dimensional ADAF, they flnd that the dynamics 
and thermal properties of the flow should be substantially changed in some parts of the parameter 
space (see also Esin 1997 for one dimensional treatment). Especially, radiative cooling and Compton 
preheating will change the polar region of the flow dramatically, and under the right conditions 
winds along the polar axis will be produced. But PO reached these conclusions based on simple 
comparison of various time scales, leaving some interesting questions unanswered. For example, 
PO showed that when Compton heating is negligible, the electrons near the polar axis cannot be 
maintained at high temperatures, suggesting the possibility of the cooling of ions. But answers to 
what temperature will electrons finally settle to and whether ions too will cool due to increased 
coupling can be obtained only after the energy equations for electrons and ions are properly solved 
with all heating, cooling, and energy exchange effects. Also, the existence of proposed solutions 
maintained by Compton preheating for the mass accretion rate at which high temperature solution 
is not possible can only be checked by solving the energy equations. 

So, in this work we explicitly integrate the ion and electron energy equations of ADAF with 
special considerations given to all relevant radiative cooling, the Compton preheating, and radiative 
transfer of Comptonizing radiation. As discussed above, we need to make fundamental assumptions 
or simplifications to deal with two dimensional nature of ADAF. Since we arc interested in the two- 
dimensional thermal properties of ADAF, we cannot use the usual height-integrated formalism. So 
we here adopt the fundamental assumption that the density and the velocity of the ADAF are 
described by the self-similar solutions of NY2, regardless of ion and electron temperature profiles. 
Since the dynamics of the flow is controlled by the pressure, and pressure, in turn, is determined 
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by energetics of electrons and ions, this approach is not completely self-consistent. In both the 
work that follows and in NY2 the pressure is dominated by ions, and since both the cooling and 
preheating that we treat affect primarily the electrons, the total pressure, which is the quantity 
most relevant to the dynamics, differs relatively little between our work and NY2. In any case, 
we believe this approach will nicely complement the extensive work on energetics of ADAF in one- 
dimensional formalism (NY3; Narayan, Kato, &; Honma 1997; Chen, Abramowicz, &; Lasota 1997; 
Nakamura et al. 1997; Gammie & Popham 1998; Popham &; Gammie 1998; Chakrabarti 1996a). 

2. Equations 

2.1. Self-Similar ADAF 

Viscous accretion flow can be heated by PdV work, viscous dissipation, and the interaction 
with radiation. It cools mostly by emission of radiation. However, locally, since the flow is moving, 
advection acts as loss of gas energy at a fixed point. The self-similar advection-dominated accretion 
flow solutions (NYl, NY2) refer to a family of solutions where viscous heating is mainly balanced by 
the advcctivc cooling with negligible radiative cooling. The dynamical property of the self-similar 
ADAF, e.g., density and velocity, is determined by the viscosity parameter a and the parameter 
e' = /~^(5/3 — 7)7(7 — 1) determined by the adiabatic index 7 of the gas and the ratio / of the 
advective cooling versus the viscous heating. When the flow is advection dominated, / = 1. 

2.2. Dynamics 

The big obstacle in investigating the axisymmetric accretion flow is the two-dimensional nature 
of the problem and the global coupling between the radiation and the gas. The physical state of the 
gas at one position is linked with all the other parts of the flow through the transfer of radiation. 
Also, the dynamics of the flow has to be determined self-consistcntly. It is very difficult at present to 
solve such two-dimensional problems in a truly self-consistent way as was done in spherical accretion 
(Shapiro 1973a,b; Park 1990a,b; Nobili et al. 1991). So we here adopt a fundamental simphfying 
assumption: the dynamics of the flow, i.e., the velocity and density of the gas, will be that of 
self-similar solutions (NY2) regardless of the thermal structure of the flow. One can alternatively 
consider this an exercise in determining the self-consistency of the NY2 solution. Thus we take 
radial velocity to be^ 

Vr{r,^)=rnk{r)v{^), (4) 



^We note, however, that recent time-dependent numerical simulations of two dimensional axisymmetric accretion 
flows very near the black hole show diverse dynamical behavior (Igumenshchev & Abramowicz 1999; Stone et al. 
1999). 
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the poloidal velocity 

= 0, (5) 

the rotational velocity 

v^ = rnk{r)n{'&)sm^, (6) 
where rQk{r) = (GM/r)~^/^. The density is 

n{r, ??)/no = m{r /rs)~^^'^n{'d), (7) 

where n^VgCTT = 1/[2(1 + ^//e)] a-i^d Thomson cross section. The free-fall spherical flow corre- 
sponds to v{d) = -v/2, ^{'d) = 0, and n(i?) = 1. As noted above, these assumptions would hold only 
if the pressure field of the final solution is similar to that of the self-similar one. We will discuss 
how self-consistent the final solutions are. 



2.3. Energy equations 

To get the relevant energy equation in axisymmetric accretion flow, we extend the relativistic 
energy equations for the spherical flow under the radiative heating and cooling to axisymmetric 
one (Park 1990a, b). Although most ADAF calculations are done in a non-relativistic framework, 
the difference in temperature between the relativistic and non-relativistic treatment is so large that 
we use here the correct relativistic energy equations. For some parameters, the output luminosity 
can be altered by more than a factor of 10. 

The ion and electron energy equations are 

Si -\- Pi dni Tms ~t~ Qie 

dr Hi dr Vr ' 

dee £e + Pe Ag - Te - Qie 



(8a) 

dr Ui dr Vr ^ ^ ^ 

where T^ig is the viscous heating rate for ions, Ag and Te the cooling and heating rates for electrons, 
£j and £e the internal energies, all per unit volume. Pi and Pe the pressures, rij and Ue the number 
densities, Tj and Tg the temperatures, of ions and electrons, respectively. This form of energy 
equations inherently accounts for the advective cooling of both ions and electrons (see Nakamura 
et al. [1997] for discussion on advection of electron energy). The energy transfer through Coulomb 
coupling is described by the rate Qie where 

3 me ,^ ,1 + 7^2(^6 + ^1 . 

Qie =7: / neniaTCikli — Kie) :== : InA, (9 

2mi^ ' ' ' ' ^ ' 072(^e + Oi)^/'' ^ ' 

Zi is ion charge, rrii and mg ion mass and electron mass, respectively, 9i = kTi/miC^, 9e = kTe/nieC^ , 
gt Thomson cross section, and In A = 20 the Coulomb logarithm. There exists a better expression 
for Qie in transrelativistic temperatures (Stepney 1983; Stepney & Guilbert 1983), but it is quite 
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costly to evaluate. The result is not sensitive to the exact functional form of Qie in the transition 
regime. 

Since electrons usually reach relativistic temperatures, = (3/2)a;*neA;Te, where x* = {2/3)[5r]6~^ — 
{r]'^-l)9-'^-l\, 7] = K;i{e-^)/K2{9-^) and K^, K2 are the modified Besscl functions (Shapiro 1973a; 
Park 1990b). The factor x* incorporates the change of the equation of state at relativistic temper- 
atures and has an asymptotic value of 1 when kTg <C mgC^ and 2 when kT^, 3> 1x1^0^. Also, ion 
energy density ei = {2>/2)x*nikTi where x* is similarly defined. The change of the equation of state 
at high temperatures can produce a significant difference compared to a nonrelativistic treatment. 
Even neglecting the derivative of x* can make an error of factor 2 (Park 1990a). Direct calculation 
of r] is numerically demanding, so we use a convenient fitting formula instead (Service 1986). 



2.4. Radiative transfer 

In a freely falling spherical accretion flow, the electron scattering optical depth from radius r 
to infinity is simply T^s(r) = m(r /rg)^^^"^ . In ADAF, the flow is moving more slowly than free-fall 
and therefore its density and optical depth is higher for the same m. We define the optical depth 
along a given polar angle 7?, 

Tesir, 7?) = m(r/r,)-V2n(^), (10) 

which has the highest value along the equator, i9 = 7r/2. For a = 0.1 and e' = 0.1, 1.0, and 10 
ADAF, n(7r/2) ~ 30, ~ 50, and ~ 400, respectively (NY2). However, n(i?) along the pole is ~ 30, 
~ 20, and ~ 10 for e' = 0.1, 1.0, and 10. Hence the flow remains optically thin along the pole for 
m < 0.1. 

It is not the intention of this work to rigorously solve the radiative transfer for the nonspherical 
flow. We add the emission from all ^ at a given r to get the averaged luminosity profile Lx{r) 
and the radiation temperature Tx{r) which is defined as 4:kTx being equal to the energy- weighted 
mean photon energy (Levich & Syunyaev 1971;Park 1990a), and use these as an approximation 
to the real radiation field. This approximation should be a good one for small e' flow where the 
density profile is almost spherical, or for low rh flow where the flow is optically thin. Moreover, 
the preheating radiation field at large radius should be well described by the streaming radiation, 
and a spherical radiation field should be good enough approximation. A more advanced treatment 
would expand the radiation field into spherical harmonics with the next (herein neglected) term 
the quadrupole component (Loeb &: Laor 1992). 

Under this simplification, Lx{t) is just the sum of all frequency-integrated emission, i.e., 
radiative cooling rate, inside radius r, 

rvout 

Lx{r) = 2TT / / Aaii{r,'&) sin M-dr^ dr. (11) 

Jrin Jo 

A far more difficult task is to calculate the Comptonization temperature Tx{r) without really 
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solving the frequency-dependent radiative transfer equation. Since we treat every emission as being 

Comptonized in situ, the main change of the radiation temperature in radius is due to the addition 
of newly produced photons to already existing photons. If the spectral energy density of existing 
photons is Ex and its radiation temperature is Tx, incremental addition of newly produced and 
locally Comptonized photons of energy density dEs and radiation temperature Ts will change the 
radiation temperature of the total to 

Ti = ^^^^=T. + dT., (12) 

where dTx is the change in Tx within dr. Solving for dTx gives (Park 1990a) 

dlnTxir) Ts{r) - Tx{r) Ks{r)r 



dlnr Tx{r) F{r) 



(13) 



where I^sif) is the i9-average of all emissions and F{r) = L{r) / A'Kr'^ . We used the approximate 
relation dEs /Ex = h-sdr/F. The T^-averaged radiation temperature, Ts{r), is calculated as 



_ Jq As {r,^)Ts{r,^) sin 



Jq As (r, 17) sm wat7 

In case where both bremsstrahlung and synchrotron emissions contribute, As{r,'t})Ts{r,'d) is re- 
placed by AcbrT^'" + Ac synT^' 



2.5. Cooling and heating 

We consider atomic cooling, bremsstrahlung, and synchrotron radiation and their inverse 
Comptonization as the main radiative cooling processes on electrons, and viscous dissipation and 
Compton scattering as the main heating processes on ions and electrons, respectively. Special care 
has been taken for absorption and inverse Comptonization of soft bremsstrahlung and synchrotron 
photons. A detailed description is given in the Appendix. 



3. Calculation 

3.1. Method of Calculation 

When the preheating radiation field is negligible, finding the solution is straightforward. One 
just integrates the energy equations (8a,b) with zero preheating luminosity. However, when the flow 
is affected by the preheating radiation, we have to find solutions that satisfy the energy equations 
(8a,b) but with a not-yet-determined radiation field. 

We use the iteration method to build the solution (Park 1990a,b). First, initial estimates for 
L°(r) = Lq and T^(r) = Tx,o are guessed for given M, m, a, e', and the ratio of the gas pressure 
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to the total (gas plus magnetic) pressure /?. The gas energy equation (8a,b) is integrated under 
this L^{r) and T^(r) for selected values of "i^'s. Now we have Ti{r,'d) and Te(r, -i?) at each {r,'d). 
Normally r is divided into 100 logarithmic intervals and •& into 10 or 20 intervals from to 7r/2. 
The angle- aver aged luminosity L^{r) and radiation temperature Tj^{r) are calculated from this 
temperature profile (Eqs. [11] and [13]). In general, Lo(ro„t) / L^{rout) and Txfi{rout) / T'Kw), 
and the solution is not self-consistent. Only for a specific combination of initial Lq and Tx-,o, will 
the model be self-consistent in the sense that the luminosity and the radiation temperature of the 
iterated solution has L}{rout) = L^{rout) and T^(rout) = T'^{rout)- We search the two-dimensional 
plane of Lq and Txfi for solutions that satisfy this self-consistency condition to within the 1 % 



There are two natural boundaries in accretion problem, inner and outer ones. At the outer 
boundary, ro^t, the state of the gas has to be specified when the gas energy equations are integrated, 
and at the inner boundary r^^, the state of the radiation field, L{rin) and Txirin), needs to be 
specified for equations (11) and (13). 

The velocity and density of gas are fixed by equations (4)- (7) for given m and M. So the only 
meaningful boundary condition is the temperature of electrons and ions. Due to the strong atomic 

cooling around 10'^ K, the gas at large radius in spherical flows (except for very low m) is most 
likely to be at the usual equilibrium temperature = Ti = T(,q K. However, self-similar 

ADAF solutions have Tg = ^ T^ir at all radii, where {5/2)kTyir = GMnip/r. So we will apply 
both boundary conditions: T{rout) = T^irifout) or Tfjcmt) = Teqiraut)- The outer boundary radius 
in most calculation is taken to be 10^7*5, which corresponds to Tyi^ ^ 10^ K. 

We adopt rin = 3rs since the radius of the marginally stable orbit around Schwarzschild black 
hole is 3r5. The luminosity at Tin would be the sum of gravitationally redshifted emission of 
radiation between rin and r^, which we approximate as 



level. 



3.2. Boundary Conditions 





(15) 



and the radiation temperature as 



Tx{rin) 




(16) 



-17- 



4. Results 

The dynamical structure of the flow, i.e., density and velocity, is fixed by the parameter e', 
which is expected to be between and 1. Esin (1997) argues that 7 = (8 — 3/3)/(6 — 3/3), and 
7 = 13/9 for P = 1/2. This would make e' = 1/2 for / = 1. But Quataert & Narayan (1999) 
argue that 7 can be close to 5/3. Also in rclativistic regime 7 approaches 4/3, and in many of 
our solutions / can be different from 1. So we choose e' = 1 as our representative value. The flow 
structure for e' = 0.1 and that for e' = 1 arc not significantly different, whereas that for e' = 10 is 
rather extreme (NY2). We also choose Yhc to be 0.1. 

4.1. Un- preheated flow 

First, we searched for solutions that have negligible preheating. This is simply achieved by 
starting the iteration with very small preheating luminosity. We first choose T{rout) = as the 
outer boundary condition, which would favor the high temperature solutions. 

For the flow onto M = lOM© or M = 10^ M© black holes, we find solutions when rh < 10~^ 
that have near virial ion temperature everywhere, which validates the self-similar ADAF. However, 
we also find that ions and electrons cool down ~ 10^ K at all radii for rh > 0.04 regardless of the 
mass of the black hole, possibly becoming cither a cool thin disk {dotted line TD in Figs. 1 & 3) or 
a cool spherical flow (triangles in Figs. 1-3). This non-existence of hot solutions agrees rather well 
with one-dimensional calculation (NY3) despite the different approach: at r = lO^r^, the critical 
mass accretion rate above which the ADAF of NY3 does not exist is m ~ 2 x 10~^ (in our units) 
for / = 0.3 and m = 10. For the intermediate mass accretion rate 10"^ ^'fn < 0.04, both ions and 
electrons near the polar axis have ~ 10^ K while the gas near the equatorial plane is maintained 
near virial temperature. As the mass accretion rate decreases, the low temperature region around 
the polar axis decreases: For the mass accretion rate m = 10~^, only the conical region of the flow 
within 6° of the polar axis is cooled, whereas for m = 10~^ that within 17° is cooled [circles in 
Fig. 5). This conflrms the fact that electrons around the polar axis will cool down, and ions too 
follow electrons due to increased Coulomb coupling between them (PO). We could not flnd any 
case where ions are kept at high temperature while electrons arc cooled down; once electrons are 
cool, the Coulomb coupling becomes too strong to allow a two-temperature plasma. These effect 
will be more severe for lower a flow because they have lower radial infall velocity. 

Figure 6 shows the typical temperature profiles of ions and electrons along different for the 
fiow with equipartition magnetic field. Solid line shows ion and electron temperature profiles for 
1? = (the polar axis) in ADAF with rh = 10"^. Up to i9 = 37r/32 (long-dashed line), the fiow stays 
near 10^ K (dotted line for ■& = 7r/32 and short-dashed line for 1? = 7r/16). Because of the very 
small viscous and adiabatic heating near the polar axis, electrons cannot be maintained at a high 
temperature, and the stronger coupling between ions and electrons at lower electron temperature 
makes ions and electrons have the same temperatures at all radii for these t^'s. Only for large 
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enough (dot-dashed lines for -& = ir/S; upper one for ion and lower one for electron) ions and 
electrons maintain their high temperatures due to higher viscous heating rate and shorter flow time. 
As result we expect a conical region with i? < 3-7r/32 around the polar axis will collapse while the 
flow along the equatorial plane accretes with significant radial velocity at high temperature. It is 
likely that the collapse would lead to an empty funnel along the polar axis. The temperature change 
across the surface of the funnel is going to be sudden, and the energy transfer by conduction may 
lead to further cooling from the hot part of the flow. But the existence and the generic shape of the 
funnel will remain the same. Of course this funnel may be filled by a tenuous outgoing wind, but 
consideration of this is beyond the scope of the present paper. However, the temperature profile 
of the flow near the equatorial plane agrees with the result from direct integration height-averaged 
equations with dynamics (Nakamura et al. 1997; Manmoto et al. 1997). Also, the flow remains at 
high temperatures for all values of mass accretion rate when the atomic cooling peak is intentionally 
removed, so it is clear that the effect we find (a collapse of the polar gas) is due to atomic cooling. 
In sum, it appears that for m > 10^^ the ADAF solutions would collapse at the pole, with the 
funnel possibly filled with a tenuous hot outgoing wind. 

We now turn to the cases where the temperature at the outer boundary is equal to the equi- 
librium temperature ~ 10^ K.^ Because of the strong atomic cooling peak near ^ 10^ K, electrons 
stay in this temperature until very strong heating takes over. For the same parameters as in 
T{rout) = Tyir cases, we find that the flow is kept at Tj = Tg ~ 10^ K for all when the mass 
accretion rate rh > 2 x 10~^ {triangles in Fig. 5). The result diflFers significantly from the prior 
boundary condition T(rout) = ^uir cases. This implies that only the very low mass accretion flow 
will reach the high temperature when the outer boundary is at ~ 10"^ K and, therefore, the cool thin 
disk would not switch to hot a ADAF automatically for intermediate mass accretion rates unless 
there is some forcing process. 

The luminosity and the radiation temperature (both at the outer boundary) of lower flow rate 
(m ^ 2 X 10~^) hot solutions for this boundary condition are shown in Figures 7 &; 8, respectively, 
as diamonds. The dimensionless luminosity / {diamonds in Fig. 7) for given dimensionless mass 
accretion rate rh is about 100 times larger than that for the spherical accretion (Shapiro 1973b) 
because the density of ADAF is roughly 20-50 times (for a = 0.1 and e' = 1), depending on {), higher 
than the spherical flow with the same rh. The radiation temperature Tx {diamonds in Fig. 8) is very 
low due to the copious soft synchrotron photons emitted (Shapiro 1973b). However, it increases as 
rh approaches 10~^'^ because now lower temperature electrons mainly cool by bremsstrahlung. The 
sudden decrease of I around m ~ 2 x 10~^ in Figure 7 reflects the transition from hot ADAF to cold 
solutions as discussed above {triangles in Fig. 5). This transition does not depend on the strength 
of the magnetic field since it is mainly determined by the atomic cooling peak near T ~ 10^ K. For 
comparison, the detailed ADAF models for Sgr A* (M = 2.5 x Mq; Narayan et al. 1998) and 
NGC 4258 (M = 3.6 x lO'^M©; Lasota et al. 1996) are shown as © and □ in Figures 7 and 8, 



^This condition would apply when (quasi) spherical flow at large radii is matched to the inner ADAF. 
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respectively. 

The luminosity and the radiation temperature of the flow without any magnetic field, thereby 
cooling only by Comptonized bremsstrahlung, are shown as stars in Figures 7 & 8. Due to the 
absence of soft synchrotron photons, the luminosity is much lower while the radiation temperature 
is much higher than the flow with magnetic field. An approximate effective temperature of the 
disk surface for the gas pressure dominated thin disk is shown as dotted line TD in Figure 8 for 
comparison. 

When conditions of self-consistency (with regard to the thermodynamics of the emitted radia- 
tion) are applied to the solutions with T{rout) = T(,q{rout) — 10^ K, we see from Figure 8 that both 
the low and high temperature branches of the ADAF solutions are only possible for rh < 10~^-^. 



4.2. Preheated flow 

PO suggested the possibility of preheated ADAF solutions, self-consistently maintained by the 
Compton preheating as in spherical accretion fiows. These solutions are found by starting from 
ad hoc initial luminosity and radiation temperature high enough to affect the flow signiflcantly 
(Park 1990a, b). Two dimensional plane of the luminosity and the radiation temperature, (I, Tx), 
is searched for the self-consistent values. 

As in un-preheated flow, the outer boundary condition is important because the physical state 
of the flow inside depends on the entropy of the flow at the outer boundary. We have two possible 
choices: T{rout) = T(,q{rout) — 10^ K or T{rout) = Tvirivout)- In most cases, the latter condition 
is awkward because T^irirout) is not the equilibrium temperature at Vout- The flow immediately 
adjusts to the equilibrium temperature Tgq 2± 10^ K and stays at that temperature until much 
stronger heating brings the temperature up. The former choice, which is more natural in usual 
spherical accretion, unfortunately, is not consistent with the self-similar ADAF, where the sum of 
the ion pressure and the electron pressure is assumed to have the virial value. However, this choice 
would make joining the ADAF solutions to outer thin disk solution much more reasonable. So we 
adopt this boundary condition. 

Figures 7 & 8 show the preheated solutions found in {m,l) and {m,Tx) plane, which may be 
called preheated ADAF (=PADAF). The circles represent the solutions with an exactly equiparti- 
tion magnetic field. We find solutions sustained by preheating in the range 6.5 x 10~^ i$ < 0.15, 
which includes the range in which no high temperature flow is possible when preheating is not 
considered (m > 0.04 or m > 2 x 10~^). The luminosity of the preheated flow increases generally 
with the mass accretion rate {circles in Fig. 7). However, for rh < 0.015 the luminosity increases 
as the mass accretion rate decreases because the electron temperature reaches above ~ 10^ K for 
this low rh flow, and synchrotron emission makes increasing contribution to the total luminosity. 
The majority of photons are now low energy synchrotron photons and the radiation temperature of 
these solutions decreases as rh decreases and I increases. This can be seen quite clearly in Figure 8. 
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The radiation temperature decrease quite rapidly as m becomes smaller than ~ 0.015. Compared 

to the (un-prehcatcd) ADAF, the PADAF has much higher radiation temperatures because most 
soft synchrotron photons are absorbed and inverse Comptonization is stronger. Solutions of NY3 
are also shown {dotted line) in Figure 7 for comparison. 

The temperature profile of the typical PADAF solution is shown in Figure 9. The parameters 
of the flow are m = 0.01, M = 10* Mq, l = A.lx 10"^, and Tx = 1.2 x 10^ K. At the outer radii, the 
flow is kept near 10^ K. When the Compton heating becomes significant, the temperature suddenly 
jumps above 10^ K first at the pole {solid line). This arises from the classic phase change due to the 
atomic cooling peak. Away from the pole, the density is higher, and the jump occurs at successively 
smaller radii. Also due to the faster infall velocity, the transition is smoother. Once the temperature 
is above > 10* K, the Coulomb coupling gets weaker, and electron temperature deviates from that 
of the ions. The flow is mostly heated by Compton preheating by hot photons produced at smaller 
radii with added help of adiabatic compression and viscous heating. The electron temperature 
increases with varying slope at Tg > 10^ K by relativistic effects and then flattens due to highly 
efficient Comptonized relativistic bremsstrahlung and Comptonized synchrotron radiation at this 
temperature range. 

Now we look at a higher m, higher I solution (Fig. 10 for rh = 0.1, M = 10 M©). Due to the 

higher density, ions and electrons are well coupled. The flow in the polar axis {solid line) is heated 
first at much larger radii compared to lower m flow, and larger i? flow is heated successively {dotted 
line for i} = ir/S, short-dashed line for i} = ir/A, long-dashed line for i} = Sir/S, and dot-dashed line 
for = 7r/2). Note that within the radius range 10^ < r/vg < 10^, the flow temperature of the polar 
region is much greater than the virial temperature {thin sohd hne) while that for the equatorial 
plane {dot-dashed line) is close to the virial value. This will produce the preheating instability 
or wind preferentially along the pole as suggested in PO because longer infall timescale along the 
pole makes this part of the flow more vulnerable to instability while shorter infall timescale near 
the equatorial plane stabilizes the flow. We may call this type of solutions ADAF with polar 
wind (WADAF). So, preheating may be yet another process to produce outflow in ADAF (Xu Sz 
Chen 1997; Blandford & Begelman 1999; Das 1999; Turolla & Dullemond 2000; see also Menou 
et al. [1999] for the polar wind in neutron star accretion). Meanwhile, if the mass accretion 
rate, and accordingly the luminosity also, is much higher, all parts of the flow will be heated well 
above the virial value. Since in spherical flow overheating outside the sonic point produces long 
term modulation on the accretion rate and that inside the sonic point produces short term flaring 
(Cowie et al. 1978), it is likely that these higher I flow may develop either a relaxational variability 
or a flaring variability. 

The range of rh to have equatorial accretion and the preheated polar outflow simultaneously 
is rather limited around m 2± 0.15. This range will be wider for larger e' and narrower for smaller 
e', because larger e' flow has a higher density contrast between the pole and the equator. 

When there is no magnetic field, Comptonized bremsstrahlung is the only cooling process, and 
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the properties of these flow show more continuous behavior {crosses in Figs. 7 &; 8). They also 
exist in a wider m range, 10~^'^ < rh < 10~^. The slope of I versus m changes around m ~ 10~^ 
because the electrons can reach near 10^ K, and relativistic radiative process becomes important. 

5. Summary and Discussion 

We have studied the two-dimensional thermal properties of the ADAF by integrating the 
energy equations under the assumption that the density and velocity of the flow is described by 
the self-similar solutions of Narayan &; Yi (1995a). Special considerations are given to the effects 
of preheating by Comptonizing hot photons produced at smaller radii and radiative transfer of 
these photons, as well as to the atomic cooling, Comptonizcd relativistic bremsstrahlung, and 
Comptonized synchrotron radiation. We find that Compton preheating is important in general. 

1. When preheating is not considered, the high temperature flows do not exist when the mass 
accretion rate is higher than ^ 10~^'^Lec~^, and even below this mass accretion rate, a roughly 
conical region around the hole cannot sustain high temperature ions and electrons. Funnels should 
exist in these flows unless rh = Mc^/Le < 10"*^. If the flow starts at large radii with a normal 
equilibrium temperature ~ 10^ K as in spherical flows, the critical mass accretion rate becomes 
m ~ 10""^-^, above which level no self-consistent solutions exist. 

2. Even above this critical mass accretion rate, the flow can be self consistently maintained at 
high temperature if Compton preheating is considered. These solutions constitute a new branch 
of solutions as in the case of spherical accretion flow; these preheated high temperature ADAFs 

can exist above the critical mass accretion rate in addition to the usual cold thin disk. Therefore, 
Compton preheating could be the mechanism to trigger the phase change from the thin disk to 
PADAF— the preheated ADAF. 

3. We also find solutions where the flow near the equatorial plane is normally accreting while 
that near the polar axis is overheated by Compton preheating, possibly becoming polar wind. The 
flow may be called WADAF — ADAF with a polar wind. 

The formalism we have used in this work is not completely satisfactory in the sense that the 
dynamics of the flow is prescribed to be ADAF like. In real flow, changes in thermal properties will 
induce those in dynamical properties. Hence, although solving the two dimensional flow structure 
with the correct radiative transfer is a daunting task, understanding the true nature of the ADAF 
like solutions seems to demand it. We feel confident, however, that the characteristic solutions for 
all accretion rates of interest (m > 10^^) will be either PADAF or WADAF type when the outer 
boundary approaches the equilibrium (T ~ 10^ K) temperature. 

We would like to thankfully acknowledge useful conversations with E. Quataert, K. Menou, 
R. Narayan, I. Yi, R. Blandford, X. Chen, and B. Paczynhski. This work was supported by Korea 
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A. Cooling and heating processes 

A.l. Atomic cooling and bremsstrahlung 

Cooling due to atomic processes and bremsstrahlung is expressed as a composite formula 
(Svensson 1982; Stepney & Guilbert 1983; Nobih et al. 1991) 

4 82 x\o-6 ) \ ' ^^^^ 

where cjt is the Thomson cross section, af the fine-structure constant, and rii the number density 
of ions. In part of the flow with high temperature where most of the radiation is produced, the 
bremsstrahlung cooling rate is not much affected by free-free absorption because most energies are 
carried by photons with energy ~ kT^ in bremsstrahlung emission. 

The cooling rate due to the Comptonized bremsstrahlung h-chr can be expressed as some 
average enhancement factor times the un-Comptonized bremsstrahlung rate 



where 



Pei = 4(^^^ 0^2(1 + 1.78101-34) foj. ^^<i 

°'-9e [ln(1.1236'e + 0.48) + 1.5] for 9e > 1 



(A3) 



27r 



Fee = -^{U-3Tr^)9^J\l + l.We + 9l-1.259l/^) fov 9e < I (A4) 
= -9e [ln(1.1230e) + 1.2746] for 9e > 1. 

TT 

Obtaining a precise estimate of the enhancement factor for arbitrary optical depth, electron temper- 
ature, and flow geometry is a quite formidable task by itself. However, Dermer, Liang, &: Canfield 
(1991) provide a convenient, yet reasonable way to deal with thermal Comptonization. They define 
a Comptonized energy enhancement factor r){i') as the average change in photon energy between 
creation and escape in a flow with an electron scattering optical depth Tgs and electron temperature 



Te, 



r]{y) = 1 + ryi - ( ^ (A5) 



where x = hiy/rUeC^, P = 1- exp(-Tes), A = 1 + A9e + lQ9j, r]i = P{A - 1)/(1 - PA), rjs = 
— 1 — In P/ In A, and 772 = rji/^^^- Although this formula is not applicable to certain ranges of u 
and Tes, it is good enough for general use. 
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We apply ri{u) to the bremsstrahlung spectrum 



€brix)d 



Abr exp - 



d 



(A6) 



to get the Comptonized bremsstrahlung cooling rate. When absorption is not important, photon 
with energy x {x < SOg) would be upscattered to r/x, and the cooling rate would be 



^Cbr = / r]X — dx+ I €br 
Jo X Jsg^ 



{x)dx 



30e 



ebr{x)dx + / {rj — l)ebr{x)dx. 
Jo 



For the second integration, we use simpler form of e^r, 



4r{x)d 



x 



X < 



X > 6e 



= Abrd 

= 

since very little emission exists above x > 9e- The resulting cooling rate is 

r/2 



A-Cbr = ^br 



1 + m- 



m + 1 



(A7) 
(A8) 

(A9) 
(AlO) 

(All) 



When 9e and t^s increases such that PA ~ 1, most photons are likely to be upscattered to 
~ fcTg. Equation (A5) has a formal divergence in this saturated Comptonization regime. If we 
assume that photons with energy above hvabs are all upscattered to kT^, the saturated energy 
enhancement factor for the whole emission is 



^Cbr/ ^br 



l+l '-^di^ 



■^abs 



X 



l + 31n 



Xabs 



Hence, the Comptonized bremsstrahlung cooling rate in general is 

V2 



K'Slr = ^br min 



l + ??2-— xT'l + 31n(— ) 



(A12) 
(A13) 

(A14) 



The energy-weighted mean photon energy for the emission is expressed as 

= Aci/ iv-r'-^dx 

Jo ^ 

- f 

Jo 



xebr{x)dx + I [r] — l)xebr{x)dx 
Jo 

1 + + 2) — h 



'^Cbr 



r?3 + 2 2r?3 + l 



(A15) 
(A16) 
(A17) 
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Since this expression has divergence near saturated regime, we put an upper hmit Tg = (Wien 
spectrum) , 



(A18) 



1 + ^?! m+i 

which has a hmit value Tg^ = when Comptonization is neghgible. 

A. 2. Synchrotron 

Even an equipartition magnetic field does not play an important role in dynamics of spher- 
ical accretion (Shapiro 1974). But synchrotron emission could be the dominant cooling process 
nevertheless. In this work we assume a magnetic field of the order of equipartition strength: 

Pm = ^= PiPg + Pm) = PP (Al9) 

where Pm is the magnetic pressure, Pg the gas pressure, and P the total pressure. 

The angle-averaged synchrotron emission by relativistic Maxwellian electrons is given by (Pa- 
cholczyk 1970) 

StT 71/ 

^syn = ^' (^Af ) dv ( A20) 

where xm = 2^'/(3i'o^g), = eS/(27rmec), and 

, 4.0505 / 0.40 0.5316\ , , o„„„ l/3^ 
I'{xm) = 1 + ^ + exp(-1.8899x;^/=^) (A21) 

is a fitting formula (Mahadevan, Narayan, &; Yi 1996). 

When absorption is not important, the cooling rate due to the optically thin synchrotron 
emission is obtained by integrating the equation (A20) 

syn 



= 213.6^ny,9l (A22) 



In generic accretion flow, a large fraction of the low energy synchrotron photons can be absorbed 
either by synchrotron self-absorption or by free-free absorption, and its effect should be incorporated 
along with that of Comptonization. 

If the flow becomes optically thick in absorption at frequency i^abs, the synchrotron cooling 
(before Comptonization) can be expressed as 



^syn — Jsyn^^syn 

where 



^syn — fsyn^syn (A23) 

j^co . poo 

fsyn= / xl'{x)dx I xl'{x)dx (A24) 
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and Xabs = huabs/fneC^- By this, we assume all photons above Vabs escape the gas flow. The 
treatment of cooling when the flow is optically thick to absorption for lower energy photons and 
thin for higher energy ones can be tricky because the absorbed radiation heats the gas and the slight 
difference in the radiation field between two different radii drives the transfer of the energy while 
higher energy photons directly escape. Correct calculation of this would require solving gas energy 
equations and radiative transfer equations with frequency dependence. A simpler way would be 
to consider only the photons that escape as the source of cooling (Ipscr & Price 1982), neglecting 
the diffusive cooling due to the absorbed photons, which we adopt here (see NY3 for different 
approach) . 

For Comptonization of absorbed synchrotron photons, we assume that photons near the cut-off 
frequency ~ I'abs are most affected by Comptonization, and the cooling rate due to the Comptoniza- 
tion of synchrotron photons is simply 

^Csyn = ri{Uabs)fsyn^%n- (A25) 

Similarly, the radiation temperature of the Comptonized synchrotron photons, Tg'"', is 

" = hUabsTjiUabs)- (A26) 



A. 3. Absorption 

There are two sources of absorption: free-free and synchrotron self-absorption. At high tem- 
perature where most of the radiation is produced, synchrotron absorption is more important in 
general. In this work, we estimate the frequency at which either absorption is important under the 
effect of Comptonization, and use the larger of the two as fobs- 

The difficulty of considering absorption under Comptonization is that a photon having a big 
chance of being absorbed in the absence of Comptonization can be upscattered and its mean 
free path for absorption can increase, resulting in escape before being absorbed. So the effective 

frequency where absorption becomes important becomes lower due to Compton scattering. This 
would be far more important to synchrotron losses than to bremsstrahlung because most of emission 
in synchrotron radiation is in the tail part of the spectrum and the precise location of the absorption 
frequency can change the emission by some exponential factor. 

Dermer, Liang, & Canfield (1991) estimate the effective absorption frequency by equating the 
characteristic timescale for absorption with that for the energy increase due to Compton scattering 
in case of free-free absorption. We apply this approach to both free-free absorption and synchrotron 
self-absorption. 

The absorption frequency u^^l^ for free-free absorption is estimated from the relation (Dermer, 
Liang, & Canfield 1991) 

VML)'^ = 1 (A27) 
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where Tff = raff and a// is the absorption coefficient 



affix) = VStt 



1 + (8/7r)V2^3/2 



■nfg 



(A28) 



and 



g = (l + 3yf/e)(l + 2^2 + 2e2)ln 



a; 



(A29) 



.2 ^.3^/2 



+(l + y4)(-^ + 2^2)^eln 
5 



47?e(11.2 + 10.4^2)^^ 



af the fine structure constant, and rg = /me(? the electron radius (Svensson 1984). This expres- 
sion is vahd for x <^ 9^ and we used y^Tr/29y^[l + (8/7r)^/^^e''^] to approximate exp(l/^e)K2(l/^e)- 
The factor T includes the effect of Comptonization, 



1 



min[l, 86e]Tt 
= 1 + ^Tes otherwise 



if y > 1 and Tgs > 1 



(A30) 
(A31) 



where ^ = 1/3 is the geometry factor and y = {AOg + 160g) max[Tes, Tg^] is the usual Compton 
parameter. 

The absorption frequency i/^^" for synchrotron self-absorption is similarly calculated (Svensson 
1984; Dermer, Liang, & Canfield 1991), 



where 



1 e^cngV I 



I'{xm) 



(Ipser & Price 1982; Mahadevan, Narayan, k Yi 1996). 

Then, the effective absorption frequency Vabs is set to be the maximum of and f^^" 



(A32) 
(A33) 



A.4. Viscous heating 

The ion heating rate per unit volume under a-viscosity description (Shakura k. Sunyaev 1973) 
is (NY2) 



Vyis = apCgflK 



3.2 (^) + p^ir})smH + {^f + sin'^i^^^^^^' 



dd 



dd 



(A34) 
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A.5. Compton heating 

We now turn to process which can be dispositive in disrupting accretion flow. Electrons in outer 
regions, can be heated by high-energy photons produced at inner hot region, 

rjrs ~ 10^ 10^, of the flow. We use Compton heating rate from the Kompaneets equation 
(Levich & Syunyaev 1971; Rybicki & Lightman 1979), 

Ak\Tx(T)-T,{T)\ Lx{t) 
1 Compt = rieCTT 2 a — T ' (Ad5j 

The information on the preheating radiation spectrum is incorporated in the angle-averaged quan- 
tities, Tx{r) and Lx{r), which are calculated from the equations (11) and (13) along with corre- 
sponding auxiliary equations. 
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Fig. 1. — Spherical and disk accretion flow solutions in (l, rh) plane [I = L/Le, rn = Mc? /Le\- 
Solid line S represents low m spherical accretion solutions of Shapiro (1973a), large crosses are those 
of Park (1990b). Low temperature solutions of high m spherical accretion are triangles (Nobili et al. 
1991) and high temperature solutions are large open circles (Park 1990b). Accretion disk solutions 
appear as dotted line TD for thin disk, dot-dashed line SD for slim disk, and dashed line NY for 
ADAF (NY3). Spherical accretion solutions with electron-positron pairs are small GUed circles 
(Park & Ostriker 1989), those with the magnetic field reconnection small squares, and those with 
preheating shocks small crosses (Chang &: Ostriker 1989). 
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Fig. 2. — The radiation temperature Tx{rout) (in K) of the self-consistent spherical solutions. 

Symbols are the same as in Figure 1. Note that the non magnetic high and low temperature 
solutions (crosses and triangles) are quite distinct on this plane, but continuous in Figure 1. The 
separate branch of the self-consistent high temperature, high- luminosity flows (circles) exists also 
in the ADAF case (PADAF). 
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Fig. 3. — Spherical and disk accretion flow solutions in {I, e) plane [e = L/Mc^]. Symbols are 
the same as in Figure 1. Note that all of the physically thick, self-consistent solutions expected or 
known to be stable {crosses, triangles, and circles) have an efficiency < 10~^. 
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Fig. 4. — Various accretion disk solutions in accretion rate m versus surface mass density S (in 
arbitrary units) plane: ADAF for ADAFs, SLE for high temperature, geometrically thin disk 
solutions, SD for slim disk solutions, RTD for radiation pressure dominated thin disk solutions, 
and GTD for gas pressure dominated thin disk solutions. 
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Fig. 5. — For the branch of the ADAF solutions for which preheating is not important, we plot the 
angle of the conical region (in degree) around the polar axis where the temperature of the flow is 
^ 10^ K as a function of tH. Civcles for boundary condition Ti^Vout) ~ T^ij. 

(rout) and triangles for 

T{rout) = Teq{rout)- With the specified cone angle pressure balance cannot be maintained and a 
conical hole would appear. Thus no solution is possible for rh above the limiting points indicated 
by the filled symbols. 
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Fig. 6. — Ion (Tj) and electron (Tg) temperature profiles of m = 0.01, M = 10^ Mq un-preheated 
flow with equipartition magnetic field for given t?'s. Soiid iine shows Ti{r;'d = 0) = Te{r;^ = 0), 
dotted line for i? = 7r/32, short-dashed hne for ■!? = 7r/16, long-dashed line for = 3-7r/32, and 
dot-dashed lines for = 7r/8 (upper one for ion and lower one for electron.) The profiles for larger 
■(J's are almost same as those of i? = vr/lG case. 
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Fig. 7. — The dimensionless luminosity L{rout)/ Le of solutions found in this work for accretion onto 
10 M0 black hole. Circles and diamonds are those with cquipartition magnetic field, and crosses 
and stars are those without magnetic field. Diamonds and stars are ADAFs where preheating is 
negligible while circles and crosses are ADAFs maintained by preheating (PADAF). Filled circle 
represent ADAF with polar wind (WADAF). One dimensional ADAF solutions (NY3) are shown 
as dotted line. Detailed ADAF models for Sgr A* (M = 2.5 x 10^ M©) and NGC 4258 (M = 
3.6 X 10'^ Mq) are shown as © and respectively. ADAFs with equipartition magnetic field are 
not permitted in the shaded region. 
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Fig. 8. — The radiation temperature Tx{rout) (in K) of solutions found in this work for accretion 
onto 10 Mq black hole with T{rout) — 10^ K. Symbols are same as in Figure 7. The effective 
temperature of the thin disk surface is shown as dotted line TD. Note that without preheating 
there are no self-consistent ADAF solutions with rh > 10~^-^. 



-40- 




Fig. 9. — Ion (Tj) and electron (Tg) temperature profiles of PADAF with equipartition magnetic 
field for given t9's: m = 0.01 and M = 10^ Mq {I = 7.6 x lO""^, Tx = 1.2 x 10^ K). Now, solid 
line shows Ti(r;'d = 0) = Te{r;'0 = 0), dotted lines for = 7r/8, short-dashed lines for = '7r/4, 
long-dashed lines for -i? = Svr/S, and dot-dashed lines for = it/2. Upper curves for ion temperature 
and lower ones for electron temperature. 
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Fig. 10.— Same as Figure 9 but for m = 0.15, M = 10 (/ = 1.0 x 10"^ Tx = 8.2 x 10*^ K) 
WADAF. Thin solid line shows the virial temperature. 



